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Abstract. The structure of the pion wave function in the relativistic constituent quark model is investigated 
in the explicitly covariant formulation of light-front dynamics. We calculate the two relativistic components 
of the pion wave function in a simple one-gluon exchange model and investigate various physical observables: 
decay constant, charge radius, electromagnetic and transition form factors. We discuss the influence of the 
full relativistic structure of the pion wave function for an overall good description of all these observables, 
including both low and high momentum scales. 

PACS. 12.39.Ki Relativistic quark model - 13.40.-f Electromagnetic processes and properties - 14. 40. Be 
Light mesons 



1 Introduction 



The understanding of the internal structure of hadrons 
within the standard model is one of the main challenge of 
nuclear and particle physics. While the high energy limit 
of the standard model will soon receive new interest from 
the expected results at LHC, a full nonperturbative de- 
scription of relativistic bound state systems in Quantum 
ChromoDynamics (QCD) is still missing. Many theoreti- 
cal frameworks already exist and shed some light on these 
systems, like QCD sum rules, lattice QCD or chiral pertur- 
bation theory. All of them have their intrinsic theoretical 
limitations. 

In order to have more physical insights into the in- 
ternal structure of hadrons, we have thus still to rely on 
constituent quark models. In the sector of up and down 
quarks, these models should be relativistic. This is also 
mandatory if one wants to understand physical observ- 
ables for which the energy scale can be large, like for in- 
stance the electromagnetic and transition form factors at 
high momentum transfer, or the decay constant of the 
pion. The interest of a phenomenological analysis of the 
structure of the pion has been renewed by recent exper- 
imental data from the Babar collaboration on the pion 
transition form factor at high momentum transfer [1] . These 
data (and older ones [2,3), as well as known data on the 
pion electromagnetic form factor [4.5.6.7.8.9. lOlfTT] and 
the precise measurement of the pion decay constant |12j 
form a rather large set of data to constrain theoretical 
models in both the low and high momentum domains. 



In the very high momentum transfer limit, factoriza- 
tion theorems enable a simple description of exclusive pro- 
cesses like the electromagnetic or transition form factors of 
the pion in terms of a distribution amplitude [13, 14, 15, 16, 
H71ITS] . This distribution amplitude is an (integrated) am- 
plitude which depends only on the longitudinal momen- 
tum fraction of the constituent quark. Corrections from 
the finite transverse momentum of the constituents may 
however contribute significantly at low and moderate val- 
ues of the momentum transfer [19120,21,22,23 . Moreover, 
the full structure of the pion involves two spin (or helicity) 
components [24] . These are a-priori of equal importance 
in this momentum range. 

The first requirement in order to build a relativistic dy- 
namical theory of bound state systems is that it should be 
invariant under the ten generators of the Poincare group. 
These generators include space time translations (four gen- 
erators), space rotations (three generators) and Lorentz 
boosts (three generators). Following this requirement, three 
forms of dynamics have been derived by Dirac already in 
1949 [25] . These are the instant form, the point form and 
the front form. 

We shall concentrate in this study on the front form. 
In this form of dynamics, the system is defined on slices 
t + = t + z — cte. This form of dynamics is of particu- 
lar interest since the boost operator along the z axis is 
purely kinematical. The electromagnetic form factors are 
thus particularly simple to calculate. However, the plane 
t + = cte is clearly not invariant under all spatial rotations. 
The angular momentum operators are therefore dynami- 
cal operators. In order to treat in a transparent way the 
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dependence of these operators on the dynamics, an explic- 
itly covariant formulation of light front dynamics (CLFD) 
has been derivedin Ref. [25]. The orientation of the light 
front plane is here characterized by an arbitrary light like 
four vector uj with uj -x = cte. This approach is a gener- 
alization of the standard light-front dynamics (LFD) [27] . 
The latter can easily be recovered with a special choice of 
the light-front orientation, uj = (1, 0, 0, —1). 

In the past few years, CLFD has been reviewed [28] 
and applied to few-body relativistic particle and nuclear 
physics. This formulation is particularly appropriate to de- 
scribe hadrons, and all observables related to them, within 
the constituent quark model. The explicit covariance of 
this formalism is realized by the invariance of the light- 
front plane uj -x — cte under any Lorentz transformation. 
This implies that uj is not the same in any reference frame, 
but varies according to Lorentz transformations, like the 
coordinate x. It is not the case in the standard formulation 
of LFD where uj is fixed to uj — (1, 0, 0, —1) in any refer- 
ence frame. Moreover, the separation of kinematical and 
dynamical transformations of the state vector provides a 
definite prescription for constructing bound and scatter- 
ing states of definite angular momentum. The dynamical 
dependence of the wave function becomes a dependence 
on the position of the light-front defined by uj. 

As we shall see in this study, this explicitly covariant 
formalism enables a very simple analysis of the structure 
of the two-body bound state. The calculation of relativis- 
tic corrections, kinematical as well as dynamical, is thus 
very easy, with a clear connection with non-relativistic 
approaches since it is also three-dimensional. A similar 
analysis in the heavy quark sector (structure of the J/'P) 
has already been done in Ref. [29] . 

In order to constrain the phenomcnological structure 
of hadron wave functions [30 ]l3Tll32] . one needs to con- 
sider several physical observables. In the case of the pion, 
this includes the decay constant, the electromagnetic form 
factor and the transition form factor. In our phenomeno- 
logical study, these observables are calculated in the rel- 
ativistic impulse approximation. Since our formalism is 
fully relativistic and can handle the full structure of the 
pion wave function - in terms of two spin amplitudes - it 
can describe low as well as high momentum scales, with 
its full kinematical structure in terms of both the longitu- 
dinal momentum fraction and the transverse momentum 
of the constituent quark and antiquark. This is at vari- 
ance with most of the previous studies which deal with 
a single distribution amplitude of the pion [T51[TH fT5, 16. 
I17|I18] which may be corrected from transverse momentum 
contributions [Til27?ll2Tll221 l2S] . A first analysis of the elec- 
tromagnetic form factor of the pion with the full structure 
of the pion wave function within the standard formulation 
of LFD can be found in Ref. [25] . 

The remainder of this paper is organized as follows. 
In section [2] we present the basic properties of CLFD. We 
apply in section[3]our formalism to the pion wave function, 
and calculate the physical observables in section [4] The 
numerical results are discussed in section[5] We summarize 
our results and present our conclusions in section [6] 



2 Covariant formulation of light-front 
dynamics 

The description of relativistic systems in CLFD has sev- 
eral nice features particularly convenient in the framework 
of the relativistic constituent quark model. The most im- 
portant ones are: 

— the formalism does not involve vacuum fluctuation con- 
tributions. Therefore, the state vector describing the 
physical bound state contains a definite number of par- 
ticles, as given by Fock state components; 

— the Fock components of the state vector satisfy a three 
dimensional equation, and the relativistic wave func- 
tion has the same interpretation as a probability am- 
plitude, like the non-relativistic one; 

— relativistic wave functions and off-shell amplitudes have 
a dependence on the orientation of the light-front plane 
which is fully parametrized in terms of the four vec- 
tor uj. In general, approximate on-shell physical ampli- 
tudes also depend on uj, whereas, exact, on-shell phys- 
ical amplitudes do not depend on the orientation of 
the light-front plane. This spurious dependence is ex- 
plicit in CLFD and is therefore under strict theoretical 
control. 

The physical bound state is described by a state vector 
expressed in terms of Fock components. The state vector 
is an irreducible representation of the Poincare group and 
is fully defined by its mass, M, its four momentum, p, 
its total angular momentum, J, and the z-axis projection 
of its angular momentum, A. The state vector, \p, A) w of 
the pion of momentum p, defined on a light-front plane 
characterized by uj (with uj ■ x = for simplicity) , is given 
in the two-body approximation by [28] 

b,A) w = (2tt) 3 / 2 £ J < os (*i,fc 2 ,p,wr) 

bl t (kr)4 2 (k 2 ) |0)*W (fcx + k 2 - p - lot) 

ol \A d3k i d3k 2 m 

2( ^ )dT (2.)3/^(2.)3/^' (1) 

where is the momentum of the quark (or antiquark) i, of 
mass m, and = -\/k 2 + m 2 . The creation operators for 
the antiquark and quark are denoted by and respec- 
tively; A is the projection of the total angular momentum 
of the system on the z axis in the rest frame and <Ji is the 
spin projections of the particle i in the corresponding rest 
system. From the delta function, 5^\ki + fc 2 — P — wr), 
ensuring momentum conservation, one gets 

V = p + ujt = k x + k 2 . (2) 

This peculiar momentum conservation law arises directly 
from the invariance of the reference system under trans- 
lations along the light-front time [28]. It is convenient to 
represent this conservation law in a systematic way. To do 
that, we shall represent on any diagram the four- vector 
ujt by a dotted line (the so-called spurion line, see [28] for 
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Fig. 1. Representation of the two-body wave function in 
CLFD. The dotted line represents the off-shell energy of the 
bound state (spurion), while the thick solid line represents the 
pion. The quark (antiquark) is shown by a thin solid line (dou- 
ble thin line). The vertex function 7~2 is defined in Eq. ( 18 1. 



more details), with an orientation opposite to the quark 
and antiquark momenta. The two-body wave function <S> 
will thus be represented by the diagram of Fig. [I] We 
emphasize that the bound state wave function is always 
an off-energy shell object (r ^ due to binding energy) 
and depends therefore on the light-front orientation. The 
parameter r is entirely determined by the on-mass shell 
condition for the individual constituents, and the conser- 
vation law ([2|. The state vector is normalized according 
to 



(3) 



The two-body wave function <&(ki, &2,p, UJT ) written 
m Eq. can be parametrized in terms of various sets 
of kinematical variables. In order to make a close connec- 
tion to the non-relativistic case, it is more convenient to 
introduce the following variables 28J defined by 



k = i- 1 ('P)k 1 =k! 



Vv 2 



n = 



uj-p 



(4) 



(5) 



where L -1 ^) is the (inverse) Lorentz boost of momen- 
tum V . The momentum k corresponds, in the frame where 
ki + k2 = 0, to the usual relative momentum between 
the two particles. The unit vector n corresponds, in this 
frame, to the spatial direction of oj. Note that this choice 
of variable does not assume that we restrict ourselves to 
this particular frame. 

The second set of variables which we shall also use 
in the following is the usual light-front set of coordinates 
(.T,RjJ, which is defined by 

x = , 

uj-p 

Rx = k% — xp , 

and where R\ is decomposed into its spatial components 
parallel and perpendicular to the direction of the light- 
front, Ri = (i? , R_l, R-ll)- We have by definition Ryto = 0, 
and thus R\ — — R^_. In the reference frame where pj^ = 0, 
Rj^ is identical to the usual transverse momentum kj_. 



The relations between these two sets of variables are given 

by 



Rl = k 2 - (n-k) 2 , 
1 _ 



.r 



(0) 



The inverse relations read: 

k 2 = R-i + r< 



n-k 



4a;(l - x) 

~R^ + m 2 

x(l — x) 



1/2 



(7) 



Note that k 2 and n-k are invariant under any rotation and 
Lorentz boost [28], like x and R^. In the non relativistic 
limit, n = n/c — > and therefore x — > 1/2 and n.k — > 0. 



3 The pion wave function 
3.1 Structure of the bound state 

The covariance of our approach allows to write down ex- 
plicitly the general spin structure of the two-body bound 
state. For a pseudoscalar particle of momentum p, com- 
posed of an antiquark and a quark of equal masses m and 
of momenta k\ and k%, respectively, it takes the form 



A=0 



1 



V2 ' 



A-, — + A* 



' uj-p 



75 v ai (kx) , (8) 



where v(ki) and u(k2) are the usual Dirac spinors, and A\ 
and A2 are the two scalar components of the pion wave 
function. For simplicity, we shall also call wave functions 
these two spin components. They depend on two scalar 
variables, which we shall choose as (cCjR^). We do not 
show for simplicity the standard isospin and color com- 
ponents of the pion wave function in Eq. ([8]). The repre- 
sentation of this wave function in terms of the variables k 
and n is given by 



1 

7=2 



fn x kl 



.91 



k 



92 W 



(9) 



where Wi are Pauli spinors and g\ t 2 are the two scalar com- 
ponents of the pion wave function in this representation. 
They depend also on two scalar variables, which we shall 
choose as (k 2 ,k.n). One can easily express A\^ in terms 
of 31,2- We get[j] 



2e,t . m 

.9i = A x H A 2 

m e k 

k A 

.92 = A 2 . 



(10) 

(11) 



1 One uses here the standard definition of 7 s with positive 
sign for its matrix elements, contrarily to where 7 s has an 
opposite sign. 
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We would like to stress that the decomposition pi is a very- 
general one for a spin zero particle composed of two spin 
1/2 constituents. In the non-relativistic limit, the com- 
ponent (?i(k 2 ,k.n) only survives and depends on a single 
scalar variable k 2 . In our phenomenological analysis, we 
shall therefore start from a non-relativistic component, 
g®, given by a simple parametrization. We shall use in the 
following either a gaussian wave function given by 



<?°(k 2 ) = a exp(-/3 k 2 ) , 
or a power-law wave function written as 

.9?(k 2 ) = 



(1 + /3 k 2 )^ ' 



(12) 



(13) 



where j3 is a parameter to be determined from experi- 
mental data, while a will be fixed from the normalization 
condition. The power 7 will be chosen equal to 2. The rel- 
ativistic component A 2l as well as dynamical relativistic 
corrections to A\ , will be calculated from radiative correc- 
tions, as explained in the next subsection. The choice (12 1 



is equivalent to the Brodsky-Huang-Lepage parametriza- 
tion [53] , 

The normalization condition writes 128 



dD <? <P 



(14) 



where dD is an invariant phase space element which can 
take the following forms, depending on the kinematical 
variables which are used 



dD 



1 



d 3 k! 



1 d 3 k 1 d 2 R_Lda; 



(2tt) 3 (1 - x)2e kl (2tt) 3 2e k (2tt) 3 2x(1 - x) ' 

With the pion wave function written in Eq. ([8]), the nor- 
malization condition writes 28 



1 = 



1 



(2tt) 3 



d 2 Rj_dx 
2x(l - x) 



Ri 



m 2 x(l — x) 



+4A 1 A 2 + 4x(l - x)A\ 



(16) 



3.2 Radiative corrections to the wave function 

In a traditional non-relativistic study of the pion wave 
function in the spirit of the constituent quark model, one 
may start directly from a simple parametri zation of the 
component g®, as given for instance in Eqs. ( 12|l3 ). How- 
ever, it is necessary to correct this wave function in some 
way in order to incorporate in a full relativistic framework 
the high momentum tail given by the one-gluon exchange 
mechanism. We shall achieve this using perturbation the- 
ory, starting from the zeroth order wave function g®. 

The (eigenvalue) equation we start from to calculate 
the bound state wave function is represented schematically 



in Fig. [2] According to the diagrammatic rules of CLFD, 
this equation writes, in the case of spin 1/2 particles [2"B] 



u(k 2 )r 2 v(k 1 ) 



d 3 ki dr' , 



8{k'i - m 2 )6>(uj-k 2 ) 



2e fcl (27r) 3 t> -ie 
xu(k 2 ) [ 7/t (# + m)r^m- #) 7l/ ] K^v(h) . (17) 

It is written in terms of the two-body vertex function F 2 
defined by [34] 



with 



u(k 2 )r 2 v(k 1 ) = (s- M 2 )$ , 



R^ + m 2 R^ + m 2 



(1-x) 



(18) 



(19) 



The mass of the pion is denoted by M w . We shall define 
for simplicity 

°-r4^ = ^( All + A2 -)^< ( 2 °) 

(s — M£) v2 V 171 U 'P J 

and similarly for O' in terms of prime quantities. The 
components A\ t2 depend on (a^R^), while A' 12 depend 

on (x,Tl' 2 ± ). 

The kernel, K^, including the appropriate color fac- 
tor, can be written as K^ v = —g^^g 2 K. in the Feynman 
gauge, with 



K = J 9 [w(fci - k[)] S [(ki - k\ + u>n - lot) 2 ] 



dri 

ri — ie 



+ / 9 [Lu{k[ - ki)} S [(k[ - ki + cjti - wt') 



where r and r' are defined by 

r= S ^ M Z, and r' = ^ 
2 oj-p 2 ui-p 

After integration over t± , we have 



T\ — ie 
(21) 



(22) 



K = 



+ 



e [ujjkx - k[)} 

-(k X - k[) 2 + 2T0J-(k 1 - k[) 
-(k[ - fcl) 2 + 2T f LU-(k' 1 - h) 



(23) 



Using the scalar products calculated in Appendix D of 
Ref. [35], one gets the following final expression for K. 

x'(x-l)e(x-x') x(x> - l)6(x> - x) 
K = F 1 F ' ( 24 ) 

where, 

/C> = m 2 (x — x')(x — x 1 — 1) 

+H A 2 x'(x - 1) + K! A 2 x(x - 1) 



(25) 



M 2 x'(x - l)(x - x') - 2x'(x - 1)Rl.R/l 
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and 

1C< 



m 2 (x — x')(x — x 
+R ± V(.t' - 1) - 
+M*x(x'-l)(x 



1) (26) 
R/l 2 .t(.t - 1) 

- x') - 2x(x' - 1)Rl.R/l . 



The quark gluon coupling constant is denoted by g, with 
g 2 = Aita s . In order to incorporate the correct short range 
properties of the quark-antiquark interaction from asymp- 
totic freedom, we shall consider in the following a running 
coupling constant a s (K 2 ), where K 2 is the off-shell mo- 
mentum squared of the gluon. It is given by K 2 = 1/JC. 
We choose a simple parametrization which gives, in the 
large K 2 limit, the known behavior given by perturbative 
QCD. We take 



a s (K 2 ) 



1 



li- 



pia»Log 



\K 2 \+A 
A 



QCD 
CD 



(27) 



At small K 2 , it is given by the parameter oP s which should 
be of the order of 1. We choose rif = 2 and Aqcd = 220 
MeV. 

In order to extract the two components Ai 2 , one should 
proceed as follows. We first multiply both sides of Eq. ( 17 1 



by u[k2) on the left and v{k\) on the right, and sum over 
polarization states. We then multiply both sides succes- 
sively by 75 and i/rf 5 , and take the trace. We end up 
with the following system of two equations for the two 
unknowns A\ 2 



Tr 



Tr 



75 (fa 
i 5 (fa + 



f m)0(fti - m) 



( S -M!)(2tt)3 
K ^d 2 R> ± dx> 



2x'(l - x 1 ) 



(28) 



and 



Tr 



Tr 



^ 5 (fa + m)0(fa-m) 
Mfa+m)A' (fa-m) 



( S -M|)(27T)3 



d 2 R^da;' 



2a/ (1 



with A' defined by 



AL 



1 ,m' 2 + m)0'{m- &) lv 



(29) 



(30) 



In perturbation theory, we shall start from a non-relativistic 
pion wave function given by 



O 



y/2 m 



75 



(31) 



where A® is calculated from ( 12 ) or (131 with g\ — g\ and 
.92 = 0. The correction to the wave function coming from 
one gluon exchange is denoted by SO and given by 





50 = -4= ( SAi — 
s/2\ m 



8A 2 



ui-p 



75 



(32) 



It is calculated from Eqs. ( 28|29 ) with the replacement 
O -> SO in the l.h.-s. and O 7 ^ O' in the r.h.-s.. The 
total wave function is then given by 

O 



o 



so , 



(33) 



with 



5Ai(x,R 2 x ) = 



SA 2 (x,Kl) = 





r d 2 R^d.T' 


(s - M 2 )2ir 2 , 


/ 2x'(l-0 


m 2 (2x' 2 -2x' 


+ 1) + R'l 


x'(l - 


x>) 


1 


f d 2 R^dx' 



(34) 



(s-M 2 )2tt 2 J 2x'{l-x') 1 V ' 
m 2 (x ~ x')(x + x' - 1) + R 2 ± x'(l - x') - R' 2 L a;(l - x) 
x(l — x)x'(l — x') 

(35) 

It is instructive to exhibit the behavior of these com- 
ponents at very high transverse momentum |Rj_|. From 
Eqs. ( 34|35 1, and using (24 1, it is easy to see that 8A\ 
is, in the absence of the running coupling constant, of the 
order of 1 /R^ while SA 2 is of the order of 1 /R^ . The run- 



R 



TP 

QCD 



ning of the coupling constant adds a factor 1 /Log 

At high transverse momentum, the relativistic component 
A 2 — 5A 2 thus dominates. More precisely, we have, in this 
limit 



A, 



AT 



3?r2 R^Log 
1 



R. 1 



QCD 



dx' 



'0 



2a;'(l-iE / ) K 



R] Log 



R 2 



Af{x) 



(36) 
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where K = x' jx for x < x' and K 
x > x' . 



(1 - x')/{l - x) for 



4 Physical observables 
4.1 Decay constant 

The pseudoscalar decay amplitude is given by the diagram 
in Fig. [3] According to the usual definition, the decay am- 
plitude is = (0\J^\tv) where is the axial current. 
Since our formulation is explicitly covariant, we can de- 
compose -T M in terms of all four-vectors available in our 
system, i.e. the incoming meson momentum p and the ar- 
bitrary position, ui, of the light-front. We have therefore 



LI 1 



(37) 



where F is the physical pion decay constant. In an exact 
calculation, B should be zero, while it is a priori non zero 
in any approximate calculation. It is a non physical, spuri- 
ous, contribution which should be extracted from the full 
amplitude 1^. Since w 2 = 0, the physical part of the pion 
decay constant can easily be obtained from 



F 



r- lj 

LU-p 



(38) 



Using the diagrammatic rules of CLFD [28], we can cal- 
culate from the graph indicated in Fig. [3] One gets, 
including color factors, 



d3kl Ar ^{kl 



m 2 )<9(wfc 2 ) 



xTr[ 



r — ie 

%7E(fc+m)r 2 (m- fa)} , (39) 



with _T 2 defined in ( 20 ) and where the notation O means 
as usual 

O = 7°0^7° . (40) 

After reduction of the scalar products, the decay constant 
is thus given by 



F 



d 2 R i dx 



2^6 

(2tt) 3 J 2x(l-x) 



[Ai + 2x(l - x)A 2 ] 



(41) 



One can immediately notice that the pion decay con- 
stant given by pi] ) is divergent with the asymptotic rel- 
ativistic component AS? given by (36). It diverges like 
Log Log Rj_/ Aq CD . This divergence is extremely soft. It 
is in fact the expression of the well known divergence of 




Fig. 3. Decay amplitude of the pion. 



radiative corrections in the process qq — > 7. To get the 
physical contribution, we just subtract the minimal con- 
tribution arising when the integral on |R_i_| is cut-off to 
Aq- The physical pion decay constant is thus 



pphys _ p _ pc 



(42) 



where 



F° 



2^6 



4 71 " 2 Log Log 




Aqcd 





dx 



2x(l 



in the limit where Aq is very large, with A^ defined in 
Eq. ([36]). 



4.2 Electromagnetic form factor 

The electromagnetic form factor is one of the most useful 
observable which can be used to probe the internal struc- 
ture of a bound state. Moreover, from the electromagnetic 
form factor at very low momentum transfer, it is possible 
to determine the charge radius of the composite particle. 
This physical observable is therefore very powerful in order 
to constrain the phenomenological structure of the wave 
function both in the low and high momentum domains. 

In the impulse approximation, the electromagnetic form 
factor is shown in Fig. [4] In CLFD, the general physical 
electromagnetic amplitude of a spinless system can be de- 
composed as [28 a 

J p = W)\e q qj p q\n(p)) = e^p+p'Y F^+^-B^Q 2 ) , 

(44) 

where e q is the charge of the quark, while is the charge 
of the pion. The physical form factor is denoted by F 7T (Q 2 ). 
In any exact calculation, B\ (Q 2 ) should be zero. We choose 
for convenience a reference frame where uj-q — 0, with 
q = p' — p. This implies automatically that the form fac- 
tors F-xiQ 2 ) and Bi(Q 2 ) depend on Q 2 — —q 2 only, since 
from homogeneity arguments their dependence on uj is of 
the form uj ■ p/ui ■ p' = 1. The physical electromagnetic 
form factor F^(Q 2 ) can be simply extracted from J p by 




Fig. 4. Pion electromagnetic form factor in the impulse ap- 
proximation. A similar contribution where the photon couples 
to the antiquark is not shown for simplicity. 
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contracting both sides of Eq. (44 1 with uj p . One thus has 

J-uj 
2 cu-p 



F^Q 2 ) 



(45) 



[28] calculated at - A 2 . We recover here naturally 
the asymptotic behavior of the pion electromagnetic form 
factor. 



By using the diagrammatic rules of CLFD, we can write 
down the electromagnetic amplitude corresponding to Fig. [4] 
where the photon interacts with the quark. Assuming it 
is pointlike, one obtains: 



F2 q {Q 2 ) 



d 3 k x dr 

2e fc , r - i 



8{kl - m 2 )0(co-k 2 ) 



dr' 



x Tr 



: S(k 2 - m 2 )9{uj-k' 2 ) 
f 2 '(k- 2 +m)-^(j£ 2 + m)r 2 (m- fa) 



2oj-p 



(46) 



where F 2 is given in Eq. ( 20 1 , and similarly for r' 2 with 
prime quantities. After calculation of the trace, one gets 



F^(Q 2 ) 



d 2 R_Lda; 
(2tt) 3 J 2x(l-x) 



+ R 2 _ -jeRx-A 

x(l — x)m 2 



2(A 1 A , 2 + A\A 2 ) + 4x{l - x)A 2 A' 2 



A X A\ 
(47) 



The wave functions A' x 2 depend on (x' , R'jJ, with x' — x 
in the impulse approximation. If we define the four mo- 
mentum transfer q by q = (go, A,qii), with A-o; = and 
qn parallel to U), we have Q 2 = —q 2 = A 2 , and thus 
R^ = R^ — iA. The contribution from the coupling of 
the photon to the antiquark can be deduced from ( 47 1 by 
the interchange x (1 — x), R± — Tl± and an 

overall change of sign. 

One thus obtains the full contribution to the electro- 
magnetic form factor of the pion 



F„(Q 2 ) = F^(Q 2 )+F™(Q 2 



(48) 



Note that this form factor, in the impulse approximation, 
is completely finite since it does not correspond to any 
radiative corrections at the jq vertex. The charge radius of 
the pion, (r 2 ) 1 / 2 , can be extracted from F n (Q 2 ) according 
to 

d 



? 2 =0 



(49) 



In the very high Q 2 limit, it is now well accepted that 
the pion form factor behaves like F ll (Q 2 ) ~ 1/Q 2 (up to 
logarithmic corrections) . This asymptotic behavior is fully 
determined by the one gluon exchange mechanism. This 
mechanism can either be considered explicitly in the hard 
scattering amplitude [13] , or incorporated in the relativis- 
tic wave function of the meson. In the spirit of the rela- 
tivistic constituent quark model, we adopt here the second 
strategy since it permits to investigate in a unique frame- 
work both low and high momentum scales. At asymptot- 
ically large Q 2 , the form factor is dominated by the con- 
tribution from the relativistic A' 2 component in Eq. ( 47 1 



4.3 Transition form factor 

The quantum numbers of the 7r transition amplitude, it — > 
7*7, are similar to the ones of the deuteron electrodisinte- 
gration amplitude near threshold, as detailed in Ref . [28^ . 
The exact physical amplitude, r p , writes therefore 



r — f p M * 



with the amplitude F pp given by 



F — -f \n v P x F 

1 pp — c^-ppvXH 1 1 17 j 



(50) 



(51) 



and where e M is the polarization vector of the final (on- 
shell) photon. The momenta P and q are defined by P = 
p + p' and q = p' — p with the kinematics indicated on 
Fig.[5j In any approximate calculation, the amplitude F^ p 
depends on uj. It should be decomposed in terms of all 
possible tensor structures compatible with the quantum 
numbers of the transition, as we did above for the decay 
constant and the electromagnetic form factor. One thus 
has HE 



F 



!!:■ — ^ £ P^fl 1 P 'F-K-f + e pt _ l y i q V L0 1 B\ + E p py-^ ' UP B 2 

+ (V p q p + V p q p )B 3 + (V p UJ p + VpU^Bi 
+ (VpPp + V pPp )B 5 , (52) 



2m 2 u)-p 



where V M = e pa/ 3 1 uj a q l3 p' 1 . From Eq. (52), we can extract 
the physical form factor by the following contraction 



F w 



2Q 2 {wp) 



Z^ vX q^xFp P 



(53) 



For the transition form factor in the impulse approxima- 
tion, the first relevant diagram, Fj} p , is indicated in Fig. [5] 




Fig. 5. Pion transition form factor in the impulse approxima- 
tion. A similar contribution where the virtual photon, denoted 
by 7*, couples to the antiquark is not shown for simplicity. 
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By applying the diagrammatic rules of CLFD, we can de- 
rive the corresponding amplitude and get 



™ _ 



2 
dr' 



{el 



4) 



d 3 kx dr 

2e 



fei 



ic 



5{k 2 - m 2 )0{uj-k2) 



le 



Tr 



- m 2 )e{uj-k' 2 ) 



(54) 



m)r 2 (m- ^i) 



The second contribution involving the coupling of the vir- 
tual photon to the antiquark can be calculated similarly. 
Other diagrams which should be taken into account at 
leading order either correspond to vacuum diagrams or 
are equal to zero for u-q = 0. After calculation of the 
trace, the total amplitude for the transition form factor 
reads 



4%/3(e 



(27T) 



d 2 R_Ldx 
2x(l - x) 



Ri 



2R ± -A + x 2 Q 2 



(55) 



A x + 2x{\ - x)A 2 - ^2~(1 - x)A 2 



The transition form factor of the pion is completely fi- 
nite thanks to the extra dependence on the transverse 
momentum as compared to the decay constant (41). The 



amplitude (54) includes a contact interaction associated 



to the elementary quark propagator between the virtual 
and real photons. It gives the factor — for' in this equa- 
tion. Additional contributions from contact interactions 
are discussed below. 



It is instructive to compare our result ( 55 ) , with the 
one obtained in the asymtotic limit using the pion dis- 
tribution amplitude [13] . This can be done by neglecting 
the mass term m 2 and the transverse momentum squared 



R_l in (55). We recover in this case the standard expres- 



sion for the transition form factor and its 1/Q 2 behavior. 
In our full calculation however, there is no need to regu- 
larize our expression in the Q 2 — > limit in order to get 
the low momentum regime |15j . 

Comparing thi s re sult with the expression of the pion 
decay constant in (41 ), one may naively identify an "equiv- 
alent" distribution amplitude given by 



d 2 R, 



1 2^/6 

F (2tt) 3 ) J 2x(l - x) 



[A x + 2x{l - x)A 2 



This "equivalent" distribution amplitude should be com- 
pared with the standard asymptotic one 4>% s = 6x(l — x) 
normalized according to 



(x)dx = 1 . 



(57) 



This however can not be done safely since the limits Q 2 — > 
oo and Rj_ 2 — > oo do not commute for the calculation 



of the transition form factor in (55). Indeed, one has to 



keep the full dependence of the transition form factor as 
a function of the transverse momentum in order to get a 
converged result in the limit Rj_ 2 — > oo. If we do the limit 
Q 2 — > oo by keeping Rj_ 2 finite, the transition form fac- 
tor is divergent, similarly to the calculation of the decay 
constant. This implies also that the "equivalent" distribu- 
tion amplitude defined in (56) is divergent when SAi and 



SA 2 are calculated from a one gluon exchange process, as 
shown in Sec. 13.21 



4.4 Contact interactions 



The contribution from one gluon exchange to the phys- 
ical observables generates in LFD several terms involv- 
ing contact interactions [2"Tll2"8] . These ones originate from 
the singular nature of the LF Hamiltonian. According to 
the diagrammatic rules given in [28], one should add, to 
each fermion (anti-fermion) propagator between two el- 
ementary vertices, a contribution of the form —fo/2u>-k 
(i/}/2uj-k) 1 where k is the momentum of the fermion. These 
contact interactions have been identified in |35j to part of 
the usual meson-exchange currents in the non-relativistic 
framework. 

For the processes under consideration in this study, we 
have thus to consider extra contributions to the pion de- 
cay constant, electromagnetic and transition form factors. 
These are shown schematically in Figs. [6| The contact in- 
teraction is indicated by a dot on these figures. One can 
easily see that since the contact interaction is proportional 
to (4 it does not contribute to the pion decay co nstant and 
electromagnetic form factor, according to Eqs. ( 38|45 ). Its 
contribution to the transition form factor is however very 
small, at most 1.5% for the highest measured momentum 
transfer. It is not included in the numerical results. 




Fig. 6. Contribution from the contact interaction (full dot) to 
the pion decay constant, electromagnetic and transition form 
factors, from top to bottom respectively, in the impulse ap- 
proximation. 
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Fig. 7. The two components of the pion wave function calculated with a gaussian parametrization in the .non-relativistic limit, 
both in the low (top curve) and high (bottom curve) momentum range. The solid line represents A\ in (31 1 while the dashed 
(dotted) line represents 8A\ (5A2) in (321. The calculation is done for x = 0.5. 




Fig. 8. Same as Fig. (|7| but for the power-law parametrization. 



5 Numerical results and discussion 

Our phenomenological analysis has three independent pa- 
rameters. The first one, /?, gives the typical size of the non- 
relativistic wave function we start from in Eqs. ( 12|13 ). 
The second parameter is the quark (or antiquark) con- 
stituent mass m. The third one is the strong coupling con- 
stant in the low momentum region given by oP s in Eq. ( 27 ) . 
The values of these parameters are indicated in Table [jj 
for the two types of non-relativistic wave functions used 
in this study. 





P 


m 


aU 


Gaussienne w.f. 


3.5 


250 MeV 


1.3 


Power-law w.f. 


3.72 


250 MeV 


0.35 



Table 1. Parameter sets of the calculation. 



These three parameters are fixed to get an overall good 
description of the pion decay constant, charge radius, elec- 
tromagnetic and transition form factors. Since the pion 



decay constant and charge radius are known with a rather 
good accuracy, we fix two of the parameters to reproduce 
these quantities (within experimental errors), while the 
third one is fixed to get an overall good account of the 
pion electromagnetic and transition form factors at mod- 
erate Q 2 . 

We show in Figs. [7]|8] the components A\ and A2 of the 
pion wave function for the two non-relativistic parame- 
trizations used in this study, in both the low and high 
momentum range. In the low momentum range, the purely 
phenomenological component A\ dominates. However, the 
contribution from one gluon exchange given by 5A\, is of 
the same order of magnitude for the gaussian parametriza- 
tion, while it is a factor 2 — 3 smaller for the power-law 
parametrization. This reflects directly the difference in the 
value of the coupling constant a®. The relativistic compo- 
nent 8A2 is always smaller, but still sizeable. 

In the high momentum domain, for Rj_ 2 > 1 — 2 GeV 2 , 
the relativistic component 5A2 dominates, as expected 
from its analytic behavior found in Sec. |3.2| We clearly see 
on these figures the interest to take into account the full 
structure of the pion wave function. It enables to describe, 



10 



O. Leitner et al.: The pion wave function in covariant light-front dynamics 



in a unique framework, both the low and high momentum 
range. 



F 


Full 


5O = 


SA 2 = 


Gaussian w.f. 


131 


92 


140 


Power-law w.f. 


131 


118 


149 



Table 2. Pion decay constant. All entries are in MeV. 



{rl) 1/2 


Full 


SO = 


SA 2 = 


Gaussian w.f. 


0.67 


0.44 


0.68 


Power- law w.f. 


0.67 


0.54 


0.68 



Table 3. Pion charge radius. All entries are in fm. 



Our predictions for the pion decay constant and charge 
radius are shown in Tables [2] and [3j respectively. The 
electromagnetic and transition from factors are shown in 
Figs. |9][l"0| for the two types of non-relativistic wave func- 
tions used in this study. Given the large experimental er- 
rors at large momentum transfer, we do not attempt in 
this study to get a best fit to all the data, but just to 
show that an overall agreement of all the available data is 
possible within our framework. 

The pion electromagnetic form factor is shown in Figs. 
[9] together with the world-wide experimental data. Given 
the experimental errors which are large above 3 GeV 2 , 
both parametrization (gaussian or power-law) give a rather 
good account of the data, in the whole kinematical domain 
available. 

In order to settle the importance of the various com- 
ponents of the pion wave function, we also show in these 
figures the electromagnetic form factor calculated with 
8A2 = (dashed line). In the kinematical domain Q 2 < 10 
GeV 2 , the contribution of the component A2 is rather 
small. This may be surprising given that A2 dominates 
the wave function for Rj_ 2 > 2 GeV 2 , as shown in Figs.[7j 
[8j This indicates that the Q 2 domain where the asymp- 
totic regime is dominant, i.e. where A2 dominates in the 
calculation of the electromagnetic form factor according 
to (47), is very high. This is in full agreement with the 
early discussions in Refs. |36j . With our numerical param- 
eters, it is above 100 GeV 2 , much above the present ex- 
perimental data. At moderate Q 2 , both low and moderate 
momentum domain of A2 dominate, and there is a partial 



cancellation between these contributions from the change 
in sign of 8A 2 at about 0.5 GeV 2 . 

The dotted line on these figures shows the contribu- 
tion of A\ only It is sizeably smaller than the full calcula- 
tion. This originates directly from the importance of the 
SAi contribution to A\ component, as shown on Figs. [7} 
[5J Note that the complete calculations using the gaus- 
sian or power-law parametrizations are extremely similar, 
the only difference being in the value of 8A\ and SA2. 



This may indicate that the most important feature, in 
this kinematical domain, is to have enough high momen- 
tum components in the pion wave function, either in the 
A\ or in theA2 components. It can come from the non- 
relativistic parametrization A\ or from the one gluon ex- 
change process giving rise to 8A\ and SA2- Since the gaus- 
sian parametrization has very little high momentum com- 
ponents, this should be compensated by a larger a®. 

The corresponding results for the pion transition form 
factor are shown on Figs.[l0] The qualitative, and to some 
extend also quantitative, features we get in this case are 
very similar to the ones detailed above for the electromag- 
netic form factor. At very high momentum transfer how- 
ever, for Q 2 > 15 GeV 2 , our results underestimate slightly 
the experimental data, with a better agreement when us- 
ing a power-law wave function. There is no way to adjust 
our parameters to get a better agreement for the transition 
form factor without spoiling the good agreement we get 
for the electromagnteic form factor. We should however 
wait for more precise experimental data before drawing 
any definite conclusions. 



6 Summary and discussion 

We have investigated in this study the full relativistic 
structure of the pion in the framework of the constituent 
quark model. This structure involves two spin compo- 
nents, which, in turn, depends on two kinematical vari- 
ables, like for instance the longitudinal momentum frac- 
tion and the square of the transverse momentum. This 
complete calculation has been made possible by the use 
of the explicitly covariant formulation of light-front dy- 
namics |28) . Our phcnomenological analysis has been com- 
pared with the full set of observables available at present: 
the pion decay constant, the charge radius, the electro- 
magnetic and the transition form factors. These observ- 
ables involve both low and high momentum scales. 

Our wave function is constructed starting from a purely 
phenomenological wave function in the non relativistic 
limit. Relativistic kinematical corrections are thus included 
exactly using CLFD, while dynamical relativistic correc- 
tions are included by a one gluon exchange process. The 
latter generates the necessary relativistic high momentum 
components in the pion wave function. 

From this full structure of the pion wave function, we 
have been able to obtain an overall very good agreement 
with all experimental data available, both in the low and 
high momentum domain. To get more physical insight into 
the relevant components of the wave function, it is however 
necessary to have more precise measurements of the pion 
electromagnetic form factor in the momentum range above 
Q 2 ~ 5 GeV 2 . It is also necessary to confirm the recent 
Babar data for the pion transition form factor at very high 
momentum transfer (till about Q 2 ~ 40 GeV 2 ), with more 
precise data. 

This analysis shows also the real flexibility of CLFD 
in describing few body systems in relativistic nuclear and 
particle physics. Its application to more fundamental cal- 
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Fig. 9. Pion electromagnetic form factor calculated with a gaussian (left plot) and a power-law (right plot) wave functions in 
the non-relativistic limit. The solid line is the complete calculation, the dotted line is the calculation without any correction 
from one gluon exchange (i.e. with 8A\ = 8A2 = 0), while the dashed line corresponds to 8A2 = 0. The experimental data are 
from P1|51[d1[71|5P1PT^PT]. 




Fig. 10. Pion transition form factor calculated with a gaussian (left plot) and a power-law (right plot) wave functions in the 
non-relativistic limit. The solid line is the complete calculation, the dotted line is the calculation without any correction from 
one gluon exchange (i.e. with 5Ai = SA2 = 0), while the dashed line corresponds to 8A2 = 0. The experimental data are from 

mm- 



culations starting from first principles is also under way 
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